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For multiqubit density operators in a suitable tensorial basis, we show that a number of nonuni- 
tary operations used in the detection and synthesis of entanglement are classifiable as reflection 
symmetries, i.e., orientation changing rotations. While one-qubit reflections correspond to antiuni- 
tary symmetries, as is known for example from the partial transposition criterion, reflections on the 
joint density of two or more qubits are not accounted for by the Wigner Theorem and are well-posed 
only for sufficiently mixed states. One example of such nonlocal reflections is the unconditional NOT 
operation on a multiparty density, i.e., an operation yelding another density and such that the sum 
of the two is the identity operator. This nonphysical operation is admissible only for sufficiently 
mixed states. 



PACS numbers: 03.65.Ud, 03.67.Mn, 03.67.-a 



The Wigner Theorem asserts that unitary and antiuni- 
tary operations exhaust all possible symmetric transfor- 
mations applicable to the wavefunction of a quantum me- 
chanical system. The unitary transformations are physi- 
cally associated with forward-in-time evolution, and an- 
tiunitary with backward-in-time evolution (see for exam- 
ple Ref. The characteristic feature of this last class 
is the presence of a conjugation operation on a wave- 
function or a transposition operation on a density op- 
erator. It is known [2| that the geometric interpreta- 
tion of the time reversing operation for a density oper- 
ator in a two-dimensional Hilbert space (aka "qubit") 
is a reflection, i.e. an orientation-changing rotation in 
0~(3) = 0(3) \ SO(3) of the corresponding Bloch vec- 
tor. 

A closely related operation, variously known as a spin 
flip , (unconditional) NOT operation, or universal in- 
verter changes the sign of the entire Bloch vector. 
In this sense it corresponds geometrically to inversion in 
the origin, which is widely known as the parity opera- 
tion For a single isolated qubit these operations are 
indistinguishable from equivalent orientation preserving 
operations, since 0(3) and SO(3) both act transitively on 
the Bloch sphere, but for multiqubit systems they corre- 
spond to partially antiunitary transformations such as 
the "partial transposition" , which can be used to detect 
bipartite entanglement 0- This highlights the intrinsi- 
cally "discrete" nature of such tests and their invariance 
under LOCC (Local Operations and Classical Communi- 
cation). 

In this paper we introduce a more general class of in- 
volutory "symmetry" operations, and argue that these 
are likewise useful in studying the multiparty nonsepa- 
rability of density operators. These operations are most 
easily described in terms of the Stokes tensor 0, IS @ 



and its "unfolding" to the so-called real density matrix 
|10|. both of which are equivalent, as carrier spaces, to 
the coherence vector [ll], [l2| . All these representations 
parametrize the real linear space of multiqubit density 
operators by the expectation values of all possible tensor 
products of the Pauli operators, differing only in their 
notations and indexing systems. The Stokes tensor in- 
dexing has the advantage of making the "affine" struc- 
ture of the set of n-qubit density operators T> n explicit, 
whereas the real density matrix has the advantages that 
both the matrix itself, as well as any operations on it 
which are diagonal with respect to the Stokes tensor, can 
be displayed as a compact 2-D table on a printed page 
(see below for examples). 

As is well-known, unitary operations on the usual Her- 
mitian density operator induces orientation-preserving 
rotations of the coherence vector, and thereby also norm- 
preserving linear group actions on the Stokes tensor / real 
density matrix [loc. cit.]. In the following, we shall fre- 
quently use the term "density" , without further qualifi- 
cation, to indicate an equivalence class of probability dis- 
tributions over an ensemble of multiqubit systems which 
all give rise to the same density operator, irrespective 
of how this is represented (as a Hermitian matrix, or a 
Stokes tensor, etc.). 

We now distinguish the following two types of nonuni- 
tary but norm-preserving operations on a multiqubit 
Stokes tensor: 

(i) local reflections applied simultaneously to two or 

more qubits; 

(ii) "nonlocal" reflections, i.e. reflections applied to the 

joint density of two or more qubits. 

The two cases are qualitatively different: while (i) is 
equivalent, up to local unitary operations, to multiqubit 
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partial transposition, (ii) is a genuinely new operation 
and does not correspond to any local operation on two 
or more qubits. In particular, the total reflection of all 
components of the Stokes tensor other than the expec- 
tation value of the identity does not correspond to time- 
reversal (i.e. to the total transpose of the density matrix) 
but rather to a multiparty NOT operation. 

Reflections on more than one qubit are nonunitary 
operations that do not necessarily yield valid (positive 
semidefinite) density operators. However, it can be 
shown that any mixed state with eigenvalues "small 
enough" is still a density operator when it is totally 
reflected. In other words, total reflection is a nonuni- 
tary involution which preserves such sufficiently mixed 
sets of density operators. On this set, total reflection 
behaves like a anti-unitary operation in the sense that 
it preserves the Hermitian structure, the trace and the 
(Hilbcrt-Schmidt) inner product. This tells us that for 
general mixed states there are more symmetries to be 
exploited than those of Wigner theorem. 

For three qubits, the set of density operators admit- 
ting a total reflection includes for example the Unex- 
tendible Product Basis (UPB) states used in Ref. 
to generate a bound entangled density operator with all 
positive partial transpositions (PPT). The "complement" 
operation that turns a separable density into the bound 
entangled UPB state is in fact a total reflection of the 
type (ii) above. The various entanglement measures (the 
concurrence, the negativity and the tangle among them) 
that rely on the use of spin-flip operations are also exam- 
ples of application of multiple one-qubit reflections of the 
type (i). In between local and total reflections lies a class 
of "nonlocal yet partial" reflections which also belong to 
the class (ii) above. These maps resemble very closely 
those used in the so-called reduction criterion [Til Il5l ] . 

Besides their unifying mathematical (group theoretic) 
character, we see reflections as a new tool to "probe" 
the structure of the set of multiparticle density opera- 
tors, in particular its nonseparable regions, by means of 
operations analogous, but inequivalcnt, to partial trans- 
position. Hopefully this will eventually lead to a bet- 
ter understanding of bound entanglement in multipartite 
systems. 



ONE QUBIT: TRANSPOSITION AND 
TIME-REVERSAL 



For a single qubit with density operator p 6 T>\ C 
C 2x2 , the Stokes tensor is the affine 3-vector [g° g T ] T , 
where g° = tr (p) /V2 = 1/V2 and g = [g 1 g 2 g 3 ] T is 
the Bloch vector of the qubit times y/2. Thus (summing 
over repeated indices) we have p — g J Xj where Aj = 
are the rescaled Pauli matrices (j — 1, 2, 3), Ao = 
and gi = tr(pXj). In this notation, the real 



density matrix is given by 



<r(p) = y/2 
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g 1 g 3 _ 
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The Stokes tensor is readily recovered from this by ap- 
plying the "col" 0| (aka reshaping 0]) operator to it 
and dividing by y/2. As is well-known, unitary transfor- 
mations of p by U G 5/7(2), namely UpU', induce rota- 
tions of the corresponding Bloch vector. This geometric 
interpretation will now be extended to antiunitary trans- 
formations [Tsl TTg| . 

Any antiunitary operation can be written as the prod- 
uct of a unitary operation and complex conjugation K. 
Given a pure state with wave vector — cq|0) + ci|l) 
(co, ci € C), let \ip) be the wave vector obtained by means 
of K alone: \^) = K\ip) = Cq\0) +c*\1). The correspond- 
ing density matrix is p = \ip)(ip\ = X) 3 -'k=o CjC%\j){k\, so 

that = EjiL3<*|j><k| = (|V>>^I) T - Since any 

density matrix is a convex combination of pure state den- 
sity matrices, the effect of K on a general p is to trans- 
pose it, i.e. p T = KpK^ = g°X + g» 1 Ai — g 2 X 2 + g 3 X 3 . 
As indicated, this is simply a change in the sign of the 
A2 component of the Bloch vector, i.e. [g 1 —g 2 g 3 ] T . 

The rotation group 0(3), of course, has two connected 
components, one of which preserves the orientation of a 
frame (namely SO(3), which contains the identity opera- 
tor I3), and one of which changes its orientation (de- 
noted here by 0~(3), to which — 13 belongs). This 
topological structure is illustrated in Fig. A reflec- 



ted) 




det(R)=-l det(R)=l 

FIG. 1: Topological structure of the rotation group 0(3). 



tion is a rotation which does not preserve orientation. 
The canonical example is spatial inversion, which is de- 
fined as multiplication by Rs = — 13- Any reflection 
R G 0~ (3) is obtained by multiplying Rs with a rota- 
tion in SO(3). For example, the reflection used in the 
transpose, Rt = diag(l, — 1, 1), can be written as the 
product of a spatial inversion with a rotation by n about 
the y-axis. 

For any vector g, spatial inversion maps g to its an- 
tipode — g = Rs{q) on a sphere of radius \\g\\ = 
\J {g 1 ) 2 + (g 2 ) 2 + (g 3 ) 2 - It follows from this together 
with the above that, for density matrices, p s = 
UKpK^W = g°X - g 1 X 1 - g 2 X 2 - g 3 X 3 where U = iX 2 G 
SU(2) rotates the Bloch vector by ir about the y-axis. In 
addition, it is easily shown that the eigenvalues of p are 
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given by 



eig(p) = / 1 



\V2\V2 



A± 11*11)} 



(2) 



Since reflections, like rotations in ^O^), are length pre- 
serving actions on the Bloch sphere, we see that the 
eigenvalues are preserved under reflections: eig(p s ) — 
eig(p T ) = eig(p). For pure states, an important differ- 
ence between R$ and Rt is that R$ maps any ket \ip) to 
an orthogonal one, whereas Rt does not. In other words, 
spatial inversion corresponds exactly to the spin-flip op- 
eration [3L l20| . 

Both the transposition and the spin-flip can also be 
defined in terms of the real density matrix, using the 
component-wise (aka Hadamard, or sometimes Schur) 
matrix product "©" . In the case of the transpose, this is 
simply: 



a(p T ) = V2 
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As shown in Ref. ^(|, an operator sum representation 
is obtained from the singular value decomposition of the 
sign matrix (left-hand factor), leading to 

a(p T ) = a(p)|0)(0|-^A 3 a(p)|l)(l| . (4) 
For the spin flip, on the other hand, it is easily seen that 



a(p s ) = V2 
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(5) 



2|0)(0|-a(p) 



These alternative representations of transposition and 
spin flip will be useful in studying multiqubit reflections 
below. 

For a single qubit the notion of reflection admits a fur- 
ther interpretation in terms of "co-completcly positive" 
(co-CP) maps. From the St0rmer-Woronowicz theorem, 
any positive 2x2 map $ is decomposable as 

$ = c$il(l-c)$ 2 oT (0<c<l), (6) 

where $1, $2 are completely positive (CP) maps and T is 
transposition. The composition $ 2 T is called a co- CP 
map. For the Bloch vector, the CP maps form a semi- 
group in the group of orientation-preserving affine maps 
GL+(3, R)©R 3 , where GL+(3, R) = {g £ GL(3, R 3 ) | 
det(g) > 0} and "(§)" denotes its semi-direct product 
with the translation group R 3 [HIU H3 Unital CP 
maps live in the GL + (3, R) component, while unital co- 
CP maps live in the other component, GL~(3, R) = {g E 
GL(3, R 3 ) I det(g) < 0}. Restricting further to sym- 
metries (i.e. trace- and norm-preserving maps), one gets 
rotations and reflections as above. 



II. TWO QUBITS: PARTIAL TRANSPOSITION, 
PARTIAL TIME REVERSAL, MULTIPLE LOCAL 
REFLECTIONS AND TOTAL REFLECTIONS 

For two qubits, a complete basis for the space of den- 
sity matrices T> 2 C C 4x4 is given by Ajk — Xj <X> Xk {j, k £ 
{0, 1, 2, 3}). This basis is also orthonormal relative to 
the Hilbert-Schmidt inner product, i.e. tr(A 3 7 £ A; m ) = 
SjiSkm for all j, k, I, m £ {0, 1, 2, 3}. For a given den- 
sity matrix p, the A-basis defines a real, rank 2 tensor Q> k 
which gives a contravariant representation of the same 
density: p = g^Ajk- Viewed as a 16- vector, g> k is affine, 



i.e. g"" = tr(pAoo) = 1/2, and it is bounded by the 15- 
dimensional sphere in R 16 of radius 1, 



tr(p 2 ) =tr(V% fe ) 



= E (^ fe ) 2 <i. 

j,k=0 



(7) 



with equality if and only if the state is pure. 

A two-qubit density matrix p is said to be separa- 
ble if it can be written as a convex combination p — 
"Ylr—l w r pA, r © PB.r for some set of real numbers w r ^ 
such that X)r=o wT = •"■> wnere PA,r> PB,r are all single- 
qubit density matrices. A necessary and sufficient con- 
dition for the separability of a two-qubit density is pro- 
vided by the positive partial transpose (PPT) criterion of 
Peres Q and Horodecki |23J. The partial transpose of a 
two-qubit density matrix p with respect to the first (left) 
subsystem A is defined as p TA = (K © I2) P (K <g> 12)^, 
and similarly p TB = (1 2 <8> K) p (1 2 ® K) . Each partial 
transpose is still a well-defined (i.e. positive semidcfinitc) 
density operator if and only if p is separable. The PPT 
criterion may be viewed as check on the feasibility of the 
"partial time reverse" operation 0, : changing the 
time arrow of one of the subsystems alone. 

In terms of the Stokes tensor g^ k , the description of 
partial transposition is very intuitive and relies on the 
observation that A2 = — X^ is the only Pauli matrix with 
imaginary elements. 

Proposition 1 For two qubits, the partial transpose op- 
erations on the density matrix p TA and p Tn act on the 
Stokes tensor (f by changing the sign of all elements 
bearing the index "2" in the corresponding subsystem: 

p TA = p ok A k + Q lk Aik- g 2k A 2 k + g 3k A 3 k (8a) 

// = O^X,,, ■ p*-A a ■ (8b) 

The verification is just a straightforward calculation, 
which may be found in Table [I] below. Note also that for 
the "total" transpose p T (— {p Ta ) Tb ) we have instead 



A jk if j, k ^ 2 or j = k = 2 (9a) 
Aj k = -A jk Hj = 2ovk = 2,j^k, (9b) 



showing that A22 behaves differently under partial or to- 
tal transposition. 



The PPT separability test of Peres-Horodecki relies 
essentially on the decomposability property ©: any 1- 
qubit positive but not CP map, when applied to a 2- 
qubit density, returns a density if and only if the original 
density is separable. Restricting from positive maps to 
symmetry operations is the same as restricting to local 
reflections. In fact, the map (|8a|l can be thought of as 
the linear transformation Rt <£> I4, where Rt is the fol- 
lowing afflne orientation-changing three-dimensional ro- 
tation: R T — diag(l, R T ) = diag(l, 1, —1,1). Since 
all single qubit reflections are unitarily equivalent, any 
matrix R £ 0~(3) can be used in place of Rt- Indeed, 
if R = diag(l, R), then local operations from the same 
connected component of 0(3) satisfy 

eig ((R ® 1 4 ) (p)) = eig ((Rt ® 1 4 ) (p)) , (10) 

where the notation must be interpreted as follows: the 
matrix R® I4 acts on the 16- vector gP k and the resulting 
16- vector provides the coefficients in the sum over the ba- 
sis elements A jk , i.e. (R®t 4 )(p) = (R <g) U) l ™p> k ^ A lm . 

Eq. I|1U|) shows that all reflections are positive but not 
completely positive. Thus we can reformulate the PPT 
criterion for the separability of two qubits as follows: 

Theorem 1 A two-qubit density matrix p is separable 
if and only if (R (g) I4) (p) is a density matrix for any 

r e o-(3). 

A particularly simple such map is Rs = diag(l, Rs), 
where Rs is the spin flip operation from the previous sec- 
tion. It is easily seen that (Rs <8> I4) (p) — 2g Qk A k — p, 
so that the sign is changed in all elements g> h except 
those appearing in the reduced density matrix of the sec- 
ond qubit (i.e. the g ok ). 

The (total) transpose p T of p corresponds to the ma- 
trix Rt i6 = Rt d> Rt = diag(l, Rt 15) with Rt 15 = 
diag(l,-l, 1, 1, 1,-1, 1,-1,-1, 1,-1, 1, 1,-1, 1)6 
50(15), where the minus signs correspond to the 6 basis 
elements obeying (|9bjl . Since the determinant of this 
matrix is positive, for two qubits the transpose is an 
orientation-preserving operation. Up to local operations 
Rt <8> Rt is equivalent to the "double local reflection" 
(or double spin flip) map Rs ® Rs- The difference be- 
tween Rs <8 I4 and Rs <8 Rs is easily understood by 
looking at Fig. |21 While Rs <S> I4 leaves the reduced 
density of the second qubit unchanged (Fig. Ufa)), the 
correlation part remains unchanged under the action of 
Rs ® Rs because its sign is flipped twice (Fig.^b)). It 
may be shown, however, that both are positive but not- 
completely-positive maps. 

All the "local" maps in 0(3) mentioned so far 
are orientation-preserving when acting on two qubits, 
even though they all have at least one factor that is 
orientation-changing when acting on a single, isolated 
qubit: det(R T ® I4) = det(R T ® Rt) = det(R s ® I4) = 
det (Rs <8> Rs) = 1- The recovery of "parity" whenever 
an orientation-changing map is applied to two or more 
qubits is due to the affine structure of the Hilbert space 
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(c) 

FIG. 2: Reflections on a 2-qubit density matrix. The two vec- 
tors contained in the smaller spheres correspond to the Bloch 
vectors g' and g ok of the two reduced density matrices, the 
third vector (double arrow) to the 2-body correlation part of 
the Stokes tensor g 3 , j, k 7^ 0: (a) the single qubit reflection 
Rs ® I4 (PPT test); (b) the double local reflection R s ® Rs 
(which is equivalent to the total transpose under LOCC); (c) 
the total reflection Rs,i6 (a nonlocal operation). 



of a qubit (resulting in a affine Bloch vector), itself a 
consequence of the trace-preserving condition: 



"1 
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'1 


0(3) 




0(3) 


50(15) 



Hence the question arises: do there exist any orientation- 
changing symmetric operations on two qubits? One such 
map is the 2-qubit total reflection Rs,m> = diag(l, — 115), 
-lis G 0(15) \ 50(15) = 0~(15). Its action (see 
Fig. Hfc)) corresponds to changing the sign of the en- 
tire tensor g^ k except for affine component g 0Q = 1/2, 
thus the name total reflection. This nonlocal operation 
is genuinely new and inequivalent to any composition of 
local symmetric operations. 

The most significant difference between total transpose 
and total reflection is that whereas the former map pre- 
serves the eigenvalues of the density matrix, the latter 
does not. Indeed, the total reflection is not even a posi- 
tive map, since it converts the density matrix of any pure 
state to one with eigenvalues [1, 1, 1, — 1]/2. This fact is 
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readily established by writing the total reflection directly 
in terms of the Hermitan density matrix as 

RsMp) = 5I4 - P, (11) 

which makes it clear that it holds for the density matrices 
of the basis states |00)(00|, |01)(01|, |10)(10| & |H)(H|, 
and that the total reflection commutes with arbitrary 
two-sided unitary transformations of p. 



The changes in the signs of the elements of the Stokes 
tensor are summarized in Table Q] for all the discrete sym- 
metric operations mentioned in this Section. It may be 
observed that Rs <8> Rs and Rt.ig — Rt ® Rt both have 
an even number of "— " signs (6), whereas Rs,w has an 
odd number (namely 15), thus confirming that a total 
reflection on a two-qubit joint density is inequivalent to 
such operations. 



TABLE I: Action (sign changes) of the rotations and reflections involving Rt, Rs and Rs,w on the components of the 2-qubit 
Stokes tensor gp . 
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# sign changes 
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12 


12 
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15 



In terms of density matrices, the positive-but-not- 
completely-positive operation ^5 Rs corresponds to 

Rs ® Rs(p) = (0-2 ® 02)9* (a 2 ® cr 2 ) = 4A 22 p*A2 2 . (12) 

The transformed density matrix p' = Rs <8> Rs (p) is fre- 
quently found in entanglement measures, including the 
concurrence C(p) = max{0, v\ — v 2 — ^3 — ^4} (where 
Vj £ eig(pp') |2J|) and the Lorentzian metric tc(pp') = 

(e 00 ) 2 - E?=i ((e 0j ) 2 + (^°) 2 ) + El k ^ k ) 2 0- ' 

III. TWO QUBITS: MATRIX STRUCTURES 
AND THE COMPUTABLE CROSS-NORM 

In this section, we show how the foregoing nonunitary 
symmetry operations on a two-qubit density matrix can 



be expressed compactly using the Hadamard product of 
matrices [2^ together with either the Stokes tensor or 
the real density matrix. We will also show that a non- 
separability criterion called the computable cross-norm 
[26j (or the matrix realignement method [2jj), which is 
inequivalent to the PPT criterion, can be computed di- 
rectly from the Stokes tensor without having to convert 
back to the traditional Hcrmitian representation. For two 
qubits, the Stokes tensor can also be viewed as a square 
array of real numbers, which is related to the real density 
matrix as follows: 
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The rearrangement of the elements seen here corresponds 
to the Choi (aka reshuffling [HI) map for n = 2 
qubits, but for n > 2 the Stokes tensor-to-real density 
matrix map is not the same as the Choi map; indeed, 
then the order of the Stokes tensor is greater than two, 
so it can no longer be identified so simply with a matrix. 

The real density matrix has the useful feature of pre- 
serving the tensor product structure of the corresponding 
Hermitian density matrix, i.e. for two qubits: a{p® p 1 ) — 
cr{p) ® cr(p') = a <£> a 1 '. It follows immediately that a 
2-qubit real density matrix can be written as a convex 
combination of 1-qubit real density matrices if and only 
if the 2-qubit density is separable. A 2 x 2 real matrix, on 
the other hand, is a real density matrix if and only if its 
upper-left element is unity and the length of the Bloch 
vector determined by the remaining elements does not ex- 
ceed unity (cf. Eq.[2J|. It should also be noted that, with 
either the real density matrix or the Stokes tensor, the 
partial trace operation involves only discarding elements 
involving the qubit traced over: no additional operations 
are needed as in the Hermitian representation. 

As shown previously for the 1-qubit case, we can 
express involutory symmetry operations by means of 
Hadamard products of the real density matrix with ma- 
trices the elements of which are all ±1. Moreover, these 
matrices will be tensor products if and only if the op- 
erations that define them are. This may be seen in the 
following list of sign matrices for all the operations given 
in Table IJ 



given below: 



R T <8> I4 
1 4 ® Rt 

Rt ® Rt 
Rs®U 
1 4 ® Rs 

Rs®Rs 



Rs. 
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(14a) 
(14b) 
(14c) 
(14d) 
(14e) 
(14f) 

(14g) 



Note that Rs,w is distinguished from the other op- 
erations not only by the fact that it is not orientation- 
preserving, but also by the fact that it is nonlocal and 
hence does not preserve the tensor product structure in 
the space of (real or Hermitian) density matrices. It is 
easily seen that the involutory mapping induced by any 
tensor product of sign matrices as above must preserve 
orientation, but there are many orientation-preserving 
mappings that are not tensor products, including the pair 



-1 +1 +1 
-1 -1 -1 
-1 -1 -1 
-1 +1 +1 



-1 +1 +1 -1 
-1 +1 -1 +1 
-1 -1 +1 +1 
-1 +1 +1 +1 



(15) 



As indicated by the double arrow, these two are related 
by the Choi map, i.e. taking the Hadamard product of 
one with the real density matrix is the same as taking the 
Hadamard product of the other with the Stokes tensor 
(cf. Eq. lO). Tests with randomly generated pure states 
quickly show that neither of these maps is positive, let 
alone completely positive. 

Similarly, one can easily construct many other discrete 
reflection symmetries which are neither locally nor uni- 
tarily equivalent to the total reflection, simply by com- 
posing the latter with any other nonlocal and nonunitary 
rotation symmetry. One interesting example is obtained 
by composing the local reflections Rs <8> Rs with the total 
reflection Rs,w on two qubits, obtaining 



{R s ®Rs)Rs,i6 



C = 



+1 +1 +1 -1 
+1 +1 -1 -1 
+1 -1 +1 -1 
-1 -1 -1 -1 



(16) 



The Hadamard product with C changes the sign of the 
bilinear (two-body) part of the Stokes tensor. It is, of 
course, a non-positive map which takes the Hermitian 
density matrix of any pure state to one with eigenvalues 
[1, 1, 1, — 1]/2. This map may also be written quite simply 
as an operator sum, as follows: 



a-^CQaip)) =■■■ 
3 

J2 {A k0 P A k0 + A ok P A ok ) - \ 1 4 • (17) 

k=l 

Finally, we show how a separability test based on the 
so-called computable cross-norm (CCN), denoted in what 
follows by , can be performed directly using the Stokes 
tensor. The CCN is a lower bound on the cross-norm 
entanglement measure in a bipartite system, denoted 
by "S", which satisfies S(p) = 1 if p is separable and 
£(p) > 1 if it is not 26]. Consequently, £(p) > 1 implies 
p is nonseparable, though not vice-versa; this condition 
is neither weaker nor stronger than the PPT criterion, 
but inequivalent to it. The CCN £ is not itself an entan- 
glement measure, since it may increase under the partial 
trace operation, but it has the advantage that it is read- 
ily computed as the sum of the singular values (aka trace 
class norm) of the reshuffled density matrix Choi(p). For 
two qubits it can also be computed directly from the 
Stokes tensor, as shown by the following: 

Proposition 2 For two qubits, the singular values of the 
Stokes tensor g ki , regarded as a matrix as in Eq. I|13|) . 
are twice those of the reshuffled density matrix Choi(p). 
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Proof. The reshuffling operation is defined to sat- 
isfy Choi(p^ (g> p 2 ) = \p2){Pi\, where \p 2 ) denotes the 
result of applying the reshaping operator to p 2 , and 
(Pi\ = \Pi) T ■ The one nonzero singular value of this ma- 
trix is simply the product of the Hilbert-Schmidt norms 
of its factors ||pi||||p2||- Recall now that p is factoriz- 
ablc if and only if the corresponding real density matrix 
a(p) is and that the linear mapping a/2 n / 2 preserves the 
Hilbert-Schmidt norm (where n is the number of qubits). 
Hence |c(p 2 )) (<r(pJ , )|/2 is the singular value decomposi- 
tion of the corresponding reshuffled real density matrix 
Choi(cr(p^ (g> p2))/2, and its nonzero singular value is 
lk(Pi)lllk(p 2 )||/2 = ||pi||||p2||. Together with the fact 
that for two qubits the Stokes tensor and the real den- 
sity matrix are related by the Choi map, this establishes 
the result for factorizable states. 

To prove the general case, we recall that the reshaping 
map Choi is self-inverse. Thus the singular value decom- 
position of a general matrix Choi(jo) = J^k £fe r fc s £ P ro ~ 
vides a canonical decomposition of p into a sum of tensor 
products J^kPkPik ® P2fc, where p k = 6= r fc II2) si |1 2 >. 
Although pk may be negative and the factors pj k , p 2 k of 
each term in this sum are not necessarily states (i.e. non- 
negative definite), we are free to apply the composition 
Choi o a to each term p\ k ® p 2 k thereby obtained. Then 
noting that er also preserves orthogonality and invok- 
ing the uniqueness of singular value decompositions com- 
pletes the proof. □ 

The claim that p is separable implies £(p) < 1 can now 
be established directly, since £(pi <g> p 2 ) = ||pi||||/02|| < 1 
and £ satisfies the triangle inequality just like any norm, 
so that for any pk > summing to unity we have 
^(J2kPkPik ® P2k) < J2kPk = 1- Tli e singular value 
decomposition of these matrices can be regarded as an 
extension of the Schmidt composition for pure states to 
mixed states. Indeed it can be shown that for pure states 
Choi(,o) has a degenerate pair of singular values which are 
equal to twice the product of the corresponding Schmidt 
coefficients. 



IV. REFLECTIONS ON THREE OR 
MORE QUBITS 

The situation is similar for three (or more) qubits, since 
the adjoint action (conjugation) still corresponds to a real 
"one-sided" rotation of the Stokes tensor, and the rota- 
tion group in all dimensions splits into orientation pre- 
serving & changing connected components. The main 
difference is that the number of inequivalent kinds of ro- 
tations and reflections goes us rapidly with the number of 
qubits. Indeed, there are 2 3 ™ local symmetries (i.e. sign 
changes in the Bloch vector components), and dividing 
this into the total number of trace-preserving discrete 
symmetries gives 

2 4"-y 2 3n = 2 4"-3n-l ^ 



locally inequivalent symmetry operations. 

It is possible, however, to identify some particularly 
significant involutions. In the case of three qubits p = 
g3 kl Ajki, the following are some of the new possibilities: 

(ia) the two-qubit partial transposition Rt,iq ® I4 (and 

the two others obtained by qubit permutation); 

(ib) the total transposition Rt,64 = diag(l, Rt, 63) 

(where i?r,63 € 50(63) is diagonal with 28 — l's 
and 35 +l's in it), which changes the sign of just 
those elements q> with an odd number of indices 
equal "2"; 

(iia) the two-qubit "reflection" Rs.m ® I4 (and the two 
others obtained by qubit permutation) - which is 
however an orientation- preserving rotation on three 
qubits; 

(iib) the total (three-qubit) reflection Rsm — 
diag(l, -1 6 3). 

The effect of Rs,64 on Q? kl is to change the sign of the 
entire tensor except for its constant component g 000 = 
l/(2y/2), showing that it may be expressed as: 

RsM = 2e 000 A oo-p - lt 8 - p. (19) 

This is again a nonlocal operation which admits no fac- 
torization into independent one-qubit operations. Simi- 
larly, the action of Rs,i6 ® I4 on g> M is to change the 
sign of the entire tensor except for the Bloch vector of 
the 1-qubit reduced density g aai . Items (ia) and (ib) 
above arc fundamentally different from (iia) and (iib). 
The first two produce a Hermitian matrix with negative 
eigenvalues whenever the density has bipartite entangle- 
ment through the cut, whereas the latter two instead 
may map even separable densities to Hermitian matrices 
with negative eigenvalues. This can be seen looking at 
the components of the UPB state. If Rs,64 is applied 
to the (separable) density puPB-scp = jEj=i l^)(V>il 

with |^) = |01+), |1 + 0), I + 01), I ) and |±) = 

(|0) ± |1)), one gets the bound entangled state pupb 
used in 01 ■ So in this case a separable state is reflected 
into an entangled state. However, no one of the 4 compo- 
nents taken alone (each is obviously separable) 
is a density when reflected. Obviously (ib) only reverses 
the time arrow on any 3-qubit density. 

In similar fashion, for any number n > 1 of qubits one 
can define an m-qubit (1 < m ^ n) nonlocal "reflection" 
Rs,4 m ® fl-4n-m, which is only a true (i.e. orientation- 
changing) reflection when the reflection is total (m = n) . 
Assuming the reflection acts on the first m qubits of an 
n-qubit density operator p, this may be written as 

R SA m®l in - m (p) = 2g°- OA »+ 1 -^Ao...oj m+ i-i» -P- 

These operations leave the norm of the n-qubit tensor 
gji---jn (1 e t r (p 2 )) invariant, but need not preserve the 
spectrum nor even leave it nonnegative, as we saw above. 
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Hence it is a "generically" ill-defined operation on the set 
of density operators of composite systems V n . 

These observations are summarized in the following: 

Proposition 3 In T> n , the linear map Rsa™ (1 < n): 

(i) preserves the trace and Hermiticity; 

(ii) preserves the Hilbert- Schmidt inner product; 

(iii) is neither unitary nor antiunitary; 

(iv) is not T> n - invariant. 

Properties (i) and (ii) together say that Rsa™ is neither 
a contraction nor a dilation map, whereas (iv) affirms 
that Rsa™ is not a positive map. 

Nevertheless, it is possible to specify a simple spectral 
condition on the density matrix that guarantees that its 
total reflection is still nonnegative definite. 

Theorem 2 Given p = g°"'°Ao...o + X <= (where \ 
is the associated homogeneous tensor), a sufficient con- 
dition for p = Rs,4 n (p) = £> "'°Ao...o — x € T^n is that the 
set of eigenvalues satisfies eig(p) C [0, 2 1 ~™]. 

Proof. The proof is based on the well-known fact [2{| 
that adding a multiple of the identity cl TO onto an m x m 
Hermitian matrix A shifts its eigenvalues by c, i.e. eig(A+ 
ct m ) = eig(A) + c. Since the eigenvalues of the random 
state's density matrix eig (g '"°Ao...o) = {2 _n } (with 
multiplicity 2"), we see that eig(p) £ [0, 2 1 ~"] implies 
both eig(x) £ [-2-", 2~ n ] and eig(- X ) £ [-2" 71 , 2"™], 
so that eig(p) G [0, 2 1 ""], as well. □ 

Hence the linear map Rs,4 n is well-defined (positive) in 
the subset of densities with eigenvalues in the interval 
[0, 2 1 -"]. 

Corollary 1 A necessary but not sufficient condition for 
Theorem\gto hold is that tr (p 2 ) < 2 1 "". 

Proof. If p Xl , . . ., p X2 „ are the eigenvalues of Xi when 
Theorem |21 holds it must be r 2 = fi 2 + . . . + p 2 ,^ < -h. 
Hence r ^ ^ and tr (p 2 ) = (g° - ) 2 + r 2 < — ^r- D 

Corollary 2 A necessary but not sufficient condition for 
p to be a density is that rank(p) 2 n ~ 1 . 

In fact, only when p is a mixture of at least 2 n ~ 1 pure 
states one can achieve eig(p) £ [0, ■^=r]- 

On a 3-qubit density, the action of Rs,ie <X> I4 is de- 
picted in Fig. |3| Essentially the entire Stokes tensor 
changes sign, except for the reduced density trAn(p)- 
Its action closely resembles the reduction criterion of 
[Til ITHj l . That criterion also makes implicit use of non- 
local reflections, but it is formulated based on a posi- 
tive map, hence it is well-posed on all of T> n . For a 3- 
qubit density it affirms that a necessary condition for 
separability is 1 2 <£> t 2 ® ^ab{p) — P ^ as well as 
I2 ® tr y i( j o) (and likewise for the other indexes). 




FIG. 3: Schematic illustration of the action of a two-qubit 
reflection -Rs,i6 ® l4(p) on three qubits. One-body (g^ 00 , 
g oko and ^ 0M ), two-body (g jk0 , g> w and g ") and three-body 
{g' kl ) correlation terms are indicated by single, double and 
triple arrows. All the signs are changed except for those of 

Q 00e . 



Since tr (I2 <8> I2 ® trAs(p) — p) = 3, one difference be- 
tween our partial reflection and the reduction criterion is 
that the latter is not a trace preserving map. Thus it is 
not a symmetry in the sense used in the paper. Nonethe- 
less, the reduction criterion utilizes a positive map which 
can be used to test separability. In our case Rs.ie ® I4 
is not a positive map even when restricted to separable 
states with eigenvalues in [0, 1/4], in which the 3-qubit 
total reflection Rs,64 is always well-posed. 

We can, however, convert our 2-qubit total reflection 
^5^6 into a "relaxed" total reflection, namely 

R r s%{p) - \{U-p), (21) 

which is the same as Rs,w applied to the "remixed" den- 
sity matrix (I4/2 + p)/3. Since the remixed density ma- 
trix now has eigenvalues in [0, 1/2], the relaxed reflection 
is a positive map by Theorem [21 It is also easily shown 
that Rg el 16 is not completely positive, and hence provides 
a necessary condition for the separability of an arbitrary 
2" x 2", n > 2, density matrix. It should be possible to 
relax all the reflections described in this paper to positive 
maps by a similar strategy. 

Concerning a total reflection, all pairs p and p = 
Rsa™{p) satisfying Theorem [21 are complementary in the 
sense that their mixture is the random state: 

\{P + P) = ^2~. (22) 

Eq. H22JI implies that Rsa™ corresponds to a multiparty 
NOT operation. In fact, also in the single qubit case, the 
NOT operation corresponds to a change of sign to the 
homogeneous part (i.e., the Bloch vector) but it is not 
modifying the sign of the trace part and hence a qubit 
and its reflection obey to J22J)- Such operation is used 
for example to map a density operator belonging to a 
subset of the Hilbert space T> n to its complement in T> n , 
for example in the UPB construction mentioned above 

mm 
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V. CONCLUDING REMARKS 

Reflections are a natural discrete class of transforma- 
tions relative to the Stokes tensor / real density matrix 
parametrization. Their meaning and relation to LOCC 
is interesting and calls for natural generalizations to non- 
local operations in the way explained above. The nonlo- 
cal reflections, in fact, originate from the nonconnected- 
ness of the group of rotations acting on the Stokes ten- 
sor parametrization. In terms of density matrices, this 
interpretation is not as sharp. As a matter of fact, op- 
erations reducible to reflections appear in the PPT test 
and in the various measures of entanglement relying upon 
"spin-flip" operations (like concurrence, negativity and 
tangle) for what concerns (multiple) 1-qubit reflections. 
Also nonlocal reflections are used: for example a total 
reflection corresponds to what is normally referred to as 
"taking the complement of a density" , used for example 
in the construction of Unextendible Product Basis states 
|l3j |. Likewise, the reduction critetion makes use of a 
positive map closely related to our nonlocal reflections. 

For the purposes of further understanding the struc- 
ture of composite quantum systems, we find it useful 
to have a unifying perspective on these nonunitary yet 
symmetric (in the sense of Wigner Theorem) transfor- 
mations. 

It is worth pointing out that reflections can be defined 
in the same terms also for SLOOC (stochastic LOCC) 
[UHH |23. For the Stokes tensor, in fact, this class of op- 



erations relaxes the group of admissible local transforma- 
tions from affine rotations to proper orthochronous local 
Lorentz transformations SO(l, 3). The reflected action in 
0~ (3) then corresponds to choosing the other connected 
component of 0(1,3) with the same time direction as 
50(1,3) (i.e., with positive "time-like" metric element). 
Also nonlocal reflections fit in with the group structure 
of nonlocal filtering operations. For example, total reflec- 
tions belong to 0(1, 4™ - 1) \S0(1, 4™ - 1). Note further 
that the idea of restricting the set of density operators in 
order to have a larger set of symmetries is "dual" to the 
idea of using group actions that are contractions [3^, . 

The idea of using reflections does not extend in a 
straightforward manner to qutrids (nor to higher dimen- 
sion quantum systems), as in this case the admissible 
parameters live on a rather complicated subset of the 
seven dimensional sphere 0, for which the rotation 
representing transposition is always admissible but spa- 
tial inversion may not be. However, the various UPB 
constructions on 3 x 3 systems of Ref. [l3j correspond to 
well-defined reflections. 

Finally notice that there are many isometries of the 
Stokes tensor that do not correspond to reflections rela- 
tive to any basis; those that do are of course involutions, 
and it would be interesting to show that any Stokes tensor 
isometry which is an involution (and hence is described 
by a symmetric orthogonal matrix) is a reflection relative 
to some basis. 
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